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A method of calculating static moments of excited states and transitions between excited states 
is formulated for non-magic nuclei within the Green function formalism. For these characteristics, 
it leads to a noticeable difference from the standard QRPA approach. Quadrupole moments of the 
first 2 + states in Sn and Pb isotopes are calculated using the self-consistent TFFS based on the 
Energy Density Functional by Fayans et al. with the set of parameters DF3-a fixed previously. A 
reasonable agreement with available experimental data is obtained. 
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I. INTRODUCTION 

To reliably predict properties of unstable nuclei at the 
modern level of microscopic nuclear theory at least two 
conditions need to be fulfilled. First, it is necessary to 
take into account the single-particle continuum which is 
especially important for the description of nuclei with 
small separation energies. Second, an approach should 
be used with the self-consistency relation between the 
mean field and effective interaction. This makes it possi- 
ble to use only one set of parameters instead of two sets, 
for the effective interaction and mean field, in non-self- 
consistent approaches. The self-consistency improved no- 
ticeably the predictive power of the theory even on the 
RPA or QRPA level, for the review see 

Nowadays, it is also necessary to add to these con- 
ditions the accounting for phonon coupling (PC). This 
problem has been studied for a long time within sev- 
eral approaches which are referred as the Quasiparticle- 
Phonon model (QPM) 0, (Q)RPA+PC 0, Extended 
Theory of Finite Fermi Systems (ETFFS) 0] (non- 
selfconsistent cases) and the self-consistent versions, 
(Q)RPA+PC [| and the ETFFS in the Quasiparticle 
Time Blocking Approximation (ETFFS(QTBA)) 0. 
For magic and semi-magic nuclei, such approaches are 
based on the fact that in these nuclei there is a small 
parameter g 2 , the dimensionless square of the phonon 
creation amplitude. For brevity, we call this weak PC 
approximation where only the g 2 terms are taken into 
account the g 2 approximation. 

In the framework of the Green function (GF) formal- 
ism, all the conditions under discussion have been re- 
alized and have shown their importance for stable nu- 




FIG. 1: g order corrections to the mass operator in magic 
nuclei. The circles with one wavy line in the first term are the 
phonon creation amplitudes. The second term is the phonon 
tadpole. 



clei too 0, H, 0j- However, all the above mentioned ap- 
proaches dealing with the PC did not take into account 
all the g 2 terms, thus limiting themselves with the pole 
diagrams only, see the first diagram in Fig. 1 where di- 
agrams for the mass operator are displayed. The second 
diagram represents the sum of all the non-pole diagrams 
that we call the phonon tadpole now. 

The problem of consistent consideration of all the g 2 
terms including tadpoles was analyzed in the article by 
Khodel [8] on the base of the general self-consistency re- 
lations for finite Fermi systems 0. The method devel- 
oped was practically realized for magic nuclei, mainly for 
static nuclear characteristics, within the self-consistent 
TFFS [10]. It was found that, as a rule, the tadpole con- 
tributions in magic nuclei are noticeable and are often of 
opposite sign as compared with those of the pole terms. 
The first attempts to include the tadpole effects for nuclei 
with pairing and for consideration of dynamical problems 
were recently made in [ill ] and [l^ . respectively. 

In the PC problem, the g 2 corrections to the mean field 
shown in Fig.l have been mainly studied up to now [lOj . 
They are, as a rule, smaller than the corresponding mean 
field values and could be partially hidden in the phe- 
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nomenological parameters used. In this work, within the 
GF method, we concentrate our attention on more del- 
icate characteristics which are proportional to g 2 them- 
selves. Namely, we analyze the static moments of excited 
states and transitions between excited states. To en- 
sure the self-consistency, we use the self-consistent TFFS 
based on the Energy Density Functional by Fayans et al. 
131 with the DF3-a set of parameters fixed previously 
14j . We briefly consider some static and low-energy char- 
acteristics in magic and generalize the method for non- 
magic even-even nuclei (Sects. II and III, respectively). 
Within this approach, and using some approximations 
we perform the first self-consistent calculations of static 
moments of the first 2 + excited states in even-even tin 
and lead isotopes (Sect. IV) 

The quadrupole moments of excited states have been 
calculated earlier within QPM in [HOil- In Ref.[l7|, the 
authors performed the calculations, which are similar to 
ours in Sect. IV, using the method later called Nuclear 
Field Theory and a phenomenological approach with the 
set of phenomenological parameters taken from experi- 
ment for each nucleus. In [l7j . a reasonable agreement 
was obtained with the experimental data for Sn and Ni 
isotopes when available at that time. The main difference 
of our approach from [Hl-Qj} is its full self-consistency on 
the (Q)RPA level and absence of any phenomenological 
or fitted parameters. These features will allow character- 
istics of unstable nuclei to be calculated. 



II. MAGIC NUCLEI 

To describe the PC effects in magic nuclei with the 
consistent account of all the g 2 terms, we follow the 
method by Khodel @ . In the g 2 approximation, the ma- 
trix element Mll for a static moment of the excited state 
(phonon) with the orbital angular moments L in a static 
external field V , is determined in terms of the change of 
the one-particle GF in the field of this phonon: 
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LL 



V°(r)5^ L G(v,r,e)dr 



de 
2m' 



where 



SlIg = S L (Gg L G) = G(e)g L G(e + u L )g L G{e) 
+G{e)g L G{e - uj L )g L G{e) + G{e)5 L g L G{e), 



(1) 



(2) 



where gh is the amplitude for the production of the L 
phonon with the energy ujl and 5l<7l, the variation of g^ 
in the field of other L phonon. Substituting Eq. (2) into 
Eq. (1), we obtain in the symbolic form: 



M LL = V°Gg L Gg L G + V°GGS L g L , 



(3) 



which is illustrated in Fig. 2 where the blocks g^ and <5l3l 
enter. 

It is convenient to transform this expression in such a 
way that the effective field V will appear instead of the 



external field V°. They are connected with the TFFS 
equation, 18] 



V = V° + FAV, 



(4) 



where T is the effective particle-hole (ph) interaction and 
A is the ph propagator (the integral over energy of the 
product of two single-particle GF's). After regrouping 
the terms in Eq.(3) and iterating the integral equation 
for the quantity 5l<7l in the effective interaction J- (for 
details, see [ToL fl2j). we obtain the ultimate expression, 



M LL = VGg L Gg L G + VGG5 L FGGg L , 



(5) 



which is illustrated in Fig. 3. It contains now the effec- 
tive field V, instead of V° in Eq. 3, and the quantity 8lF 
in the second term which denotes the variation of the ef- 
fective ph interaction T in the field of the L phonon. For 
the density dependent TFFS effective interaction F{p), 
the following ansatz can be readily obtained [1, [l(| : 



3T 

8lF(t) = -fpt(r)Y LM (n), 



(6) 



where p l £ = Ag^ is the transition density for the L 
phonon excitation. The first term of Eq.(5) coincides 
with the result by Speth [T|| while the second one with 
the 6l T quantity is a generalization to take into account 
all the g 2 terms. 



M LL 




FIG. 2: Matrix element Mll for magic nuclei, Eq. 




FIG. 3: Matrix element M LL in the form of Eq. J5j. 

All the above equations can be readily modified for 
such processes as the transition between the excited 
states L and V in the external field V°(uj = ujl' — <^l) 
or the excitation of the two-phonon state L + L' in the 



external field V (oj — uj^i 
case corresponds to u> = 0. 



ujl)- The static moments 



III. NON-MAGIC NUCLEI. COMPARISON 
WITH QRPA. 

In the case of nuclei with pairing, it is necessary to use 
four GF's (G,G h ,FW,FW m the usual notation [3). 
To describe phonons, one has to use the complete set of 
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the QRPA equations which include the ph, hp, pp, and 
hh channels and four effective fields V, V h . d (1 ) and rf' 2 ' 
[l8j , respectively. As the pp and hh channels give a small 
contribution in the case of the first 2 + levels [2fJ , which is 
considered in the next Section, we do not consider these 
channels and, accordingly, the fields S 1 ^ and S- 2 \ Then 
we obtain eight terms M LL , instead of one in Eq.(5). The 
typical two terms , and are shown in Fig. 4. 

In this study we will consider the case without the terms 
with SlJ 7 and SlJ 7 ^, see Sect. IV 

The integral of the three GF's A^2 3 a " has the form 
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J l v 2 "3 
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2 2 2 



n u 2 "3 



Ul + E\2 \E 23 — UlL' Ei 3 +UJl 
1 



(7) 



j x u, 2 v 3 



J l"2 Lt 3 



UJ — E\2 \£'23+ w L' Ei 3 — UJl 



where E 12 = E x + E 2 ,E 1 = y /(e 1 - fi) 2 + A\ and the 
low index 1 = (ni,li,j\) (spherical nuclei) is the set of 
single-particle quantum numbers. 

In the pairing case there are eight such integrals of 
three GF's ^i23 air , where i = 1-8. After a long algebra, 
one can obtain the final formula for the diagonal matrix 
element Mll with L = L', which corresponds to the case 
of static quadrupole moment of the excited L state. Just 
this case is considered in the next section . 
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Let us compare this expression with the respective re- 
sult of QRPA approach. Here we mean the usual way 
which uses the QRPA wave functions for the matrix el- 
ement between two excited states. In Ref.[2l| the ex- 
pression for the matrix element (the B(E2) quantity, to 
be exact) has been derived using the bare external field 
and, like in our case, the QRPA wave functions without 
the pp and hh-channels. The first square brackets (in the 
first half of Eq. (9)) coincide completely with the factor 
Vi2 u 23 u ti m Refs.0, Thus, the first half of Eq.(9) 
corresponds to the expression v^ 2 (V'23'03i+ ( ^23</ ) 3i) in [21| 
because the phonon amplitudes ip and <fi contain, by def- 
inition, the denominators (E\ 2 ~ojl) and (E\ 2 + re- 
spectively. Therefore, the second half of Eq. (9) (with the 
common factor \jE\ 2 \ generalizes the usual QRPA ap- 
proach. This part of Eq.(9) describes the contribution of 
the ground state correlations (the so-called graphs going 
back) to our "triangle" with the integral of three GF's, 
see Eq. (7). We will calculate the quantitative contribu- 
tion of such correlations in the next section. The second 
generalization is the appearance of the effective field V, 
which depends on the frequency uj = lol ± w' L instead 
of external field V°, which does not depend on the fre- 
quency. The terms with 5lF and 8lF^ are the third 
generalization of the QRPA approach 



IV. CALCULATIONS OF STATIC 
QUADRUPOLE MOMENTS OF THE FIRST 2+ 
STATES IN TIN AND LEAD ISOTOPES 

The quadrupole moment of the excited state L is con- 
nected with the matrix element Mll, Eq. (8) (I = 
2, V(r) = V~(r)F 20 (n) ), as follows: 



Q 



16tt 



M 



LL- 



(10) 



In the recent work [22|], we calculated the static 
quadrupole moment of the first 3~ level in 208 Pb, using 
only the first term in Eq. (5) (or Fig. 3), and obtained 
the value Qtheor = ~~ 0.26 b which is in a reasonable agree- 



ment with the experimental one, Q c 



0.34±0.15 b 



[23l |. At present, we have calculated the second term 
with SlJ 7 from Eq. (6) and found that it gives approxi- 
mately 1% as compared with the first term. Though the 
contribution of the terms with 8lF and 8lF^ in non- 
magic nuclei should be analyzed specially, it can hardly 
be considerable in the problem under consideration and 
is omitted in the present calculations. 

We calculated the quadrupole moments of the first 2 + 
states in non-magic tin and lead isotopes according to 
Eqs. (8 - 10) in the A-representation with self-consistent 
single-particle wave functions tf>\ obtained within the 
EDF method [TJ] with the functional DF3-a [14]. The 
spherical box of the radius i?=16 fm was used to sim- 
ulate the single-particle continuum. We examined the 
dependence of calculation results on the cut-off energy 




ensures 1% accuracy. To calculate t he q uantities V and 
gL, the results of our previous article [20( have been used 
where all the calculations were performed in the coordi- 
nate representation using the same self-consistent DF3-a 
basis as in the present calculation of the matrix element 
Mll- Thus, no fitted parameters was used in the present 
calculations. 

The results are given in Table 1 and Figs. 1, 2. Ex- 
cept for the case of 112 Sn and 208 Pb nuclei, we obtained 
a reasonable agreement with experimental data [23| . Un- 
fortunately, they have rather big errors so that we need 
some better measurements to check our approach. 

We have also calculated the contribution of the ground 
state correlations term in Eq. (9) (see the discussion 
at the end of Sect. Ill) and obtained that it is rather 
considerable. For some nuclei it is about 60% of all the 
triangle contribution. This very interesting result will be 
discussed in more detail separately. 
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FIG. 5: (Color online) Quadrupole moments of the first 2+ 
excited states in even Sn isotopes. 



TABLE I: Quadrupole moments Q (e b) of the first 2 + states 
in Sn and Pb isotopes. (Q n and Q p are the neutron and 
proton contributions to the final result Qtot~Q n + Q v ■ ) 



nucl. 


Q n 




Qtot 




lui Sn 


-0.05 


-0.01 


-0.07 


— 


104 o 

Sn 


-0.18 


-0.03 


-0.21 




106 Sn 


-0.28 


-0.05 


-0.33 




108 Sn 


-0.31 


-0.07 


-0.38 




110 Sn 


-0.38 


-0.10 


-0.48 




112 Sn 


-0.32 


-0.11 


-0.43 


-0.03(11) 


114 Sn 


-0.15 


-0.11 


-0.26 


0.32(3), 0.36(4) 


116 Sn 


0.00 


-0.10 


-0.10 


-0.17(4), +0.08(8) 


118 Sn 


0.10 


-0.09 


0.01 


-0.05(14) 


120 Sn 


0.12 


-0.08 


0.05 


+0.022(10), -0.05(10) 


122 Sn 


0.09 


-0.07 


0.02 


-0.28 < Q <+0.14 


124 Sn 


0.00 


-0.06 


-0.06 


0.0(2) 


126 Sn 


-0.08 


-0.05 


-0.12 




128 Sn 


-0.10 


-0.03 


-0.13 




130 Sn 


-0.05 


-0.01 


-0.07 




132 Sn 


0.04 


0.00 


0.04 




134 Sn 


0.00 


-0.01 


-0.01 




190p b 


-0.60 


-0.29 


-0.89 




192p b 


-0.77 


-0.35 


-1.12 




194p b 


-0.90 


-0.39 


-1.28 




196 Pb 


-0.85 


-0.38 


-1.23 




198 Pb 


-0.67 


-0.35 


-1.02 




200p b 


-0.27 


-0.23 


-0.50 




202 Pb 


0.02 


-0.15 


-0.13 




204p b 


0.18 


-0.07 


0.11 


+0.23(9) 


206 Pb 


0.11 


-0.02 


0.10 


+0.05(9) 


208 p b 


0.01 


0.04 


0.05 


-0.7(3) 



V. CONCLUSION 

We have considered the method to calculate static 
moments of the excited states and transitions between 
excited states, which, generally speaking, are described 
within the QRPA, with taking all the g 2 terms into ac- 
count in magic and non magic nuclei. It was shown that, 
in addition to the old results [13, Qjl, new terms with 



5 




-1,4 I — ' — I — ' — I — ' — I — ' — I — ' — I — ' — I — ' — I — ' — I — ' — I — ' — I — ' — I 

188 190 192 194 196 198 200 202 204 206 208 210 

A 



FIG. 6: (Color online) Same as in Fig. 5 but for Pb isotopes 

5t,F and S^F^ appear, which contain the density de- 
pendence of both the ph and pp effective interactions. 
We have performed the self-consistent calculations of the 
static quadrupole moments of the first 2 + states for Sn 



and Pb isotopes using the known functional's parame- 
ters set DF3-a. Except for the 112 Sn and 208 Pb cases, 
a reasonable agreement has been obtained with the ex- 
periment available, though, as a rule, the experimental 
data have big errors. In these calculations we did not 
take into account new terms with 6lF and SlF^ . As 
all modern microscopic calculations deal with the den- 
sity dependence of ph and pp effective interactions, it is 
important to analyze their role in the problems under 
consideration. This rather complicated problem will be 
considered separately. 
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